Non-universal BCS-BEC crossover in resonantly interacting Fermi gases 
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We investigate the non-universal behavior of the BCS-BEC crossover model at the normal to su- 
perfluid transition. By using a finite temperature quantum field theoretical approach due to Nozieres 
and Schmitt-Rink and by making the effective range expansion of the effective two-body interac- 
tion we numerically calculate the crossover transition temperature as a function of the scattering 
length, ap, and the effective range parameter, r^. In an ultracold two-component atomic Fermi gas 
BCS-BEC crossover physics is expected to appear near magnetic-field-induced Feshbach resonances. 
By matching two-body scattering properties near a Feshbach resonance to a simple renormalized 
model potential, a broad resonance is characterized by a small effective range and the gas displays a 
universal BCS-BEC crossover behavior. On the other hand, for a narrow resonance thermodynamic 
quantities depend the effective range which may be large and negative at resonance. For increas- 
ing values of — rc we find the transition temperature to be suppressed in the crossover region. We 
furthermore argue for existence of a lower bound of Tc at fixed coupling for increasing — rc in the 
crossover region. 

PACS numbers: 03.75.Ss, 03.75.Hh, 05.30.Jp 



I. INTRODUCTION 

The prospect of creating a fermion superfluid in 
trapped atomic Fermi gases using a magnetic field in- 
duced Feshbach resonance has received much experimen- 
tal and theoretical attention in the past year. Recently 
long-lived diatomic molecules of fermionic ^Li and ""^K 
Bose-Einstein condensation (BEC) of molecular 
^Li2 and '*°K2 dimers crossover from a molecular 

BEC to a degenerate and strongly interacting Fermi gas 
0,01, and condensation of fermionic atom pairs [fll 
have been reported. From measurements of the energy 
gaps in the excitation spectrum 12] and of the fre- 
quencies and damping of compressional collective modes 
as a function of the effective atom-atom interaction [3.11^ 
evidence is emerging that some of the experiments are al- 
ready probing the superfluid phase. 

From a theoretical point of view the low temperature 
atom gases with resonant interactions have on the one 
hand been described in terms of coupled boson-fermion 
models where the atomic and molecular degrees of free- 
dom have been treated separately but constrained by to- 
tal fermion number conservation. On the other hand, the 
description of superfluid pairing in Fermi systems with a 
varying coupling constant has also been analyzed using a 
generalized Bardeen-Cooper-Schrieffer (BCS) variational 
function at zero temperature and b y G reen's func- 
tion methods at nonzero temperatures [ifll ll^ includ- 
ing explicitly only the fermionic degrees of freedom. In 
the latter approach a boson-like two-body bound state 
is formed at large coupHngs and the system displays a 
smooth transition from a gas of fermions which under- 
goes a phase transition to a BCS superfluid to a gas of 
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bosons which condenses at the Bose-Einstein condensa- 
tion temperature. 

In ^3,0,01 it was argued that the universal crossover 
model based on a static contact potential appHes only 
to broad Feshbach resonances, and that for more gen- 
eral potentials {e.g., potentials characterized by the pres- 
ence of a repulsive barrier) the solutions to the crossover 
equations depend sensitively on the effective range, which 
gives a measure of the energy dependence of the scatter- 
ing phase shift and transition matrix. Similar conclusions 
were drawn in ^2Q| based on a microscopic flve coupled 
channels calculation of diatomic ^Li. In this study the 
molecular Born-Oppenheimer potentials were adjusted 
in order to reproduce the experimental resonance posi- 
tions of the broad resonance at 822 G, the narrow res- 
onance at 543G in addition to the zero crossing of the 
scattering length at 533 G and the hyperflne splitting 
between = 3/2 and F — 1/2 states. The effective 
range was extracted and found to be small and posi- 
tive for the broad resonance and thus consistent with the 
standard BCS-BEC crossover model. In contrast, the ef- 
fective range for the narrow resonance was found to be 
large and negative. In relation to the two-channel scat- 
tering model derived from the boson-fermion models the 
relation between the effective range, at resonance and 
the boson-fermion coupling constant, gsF, was found to 
be, To — — 87r?i'*/(m^|5s_Fp), where m is the mass of the 
atoms and when neglecting the background scattering 
length Q. 

In the present work we generalize the calculation of the 
BCS-BEC crossover transition temperature to potential 
models with nonzero effective range. The paper is orga- 
nized as follows: Section UTI contains a description of the 
flnite temperature Green's function formalism [l5| due 
to Nozieres and Schmitt-Rink (NS). Partly following 
the thermodynamic equations determining the transition 
temperature and chemical potential are presented, and 
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the renormalization of the T-matrix equation including 
the effective range expansion of the effective interaction 
is explained. In Section IIIII the resulting equations are 
numerically solved to determine the critical chemical po- 
tential and transition temperature in the crossover region 
as a functions of the effective range. 



II. BCS-BEC CROSSOVER MODEL 

The BCS-BEC crossover problem has a long history 
0,^5, 16, 22] and it was recently revived in the context 
of high-Tc superconductivity of the cuprates. The generic 
crossover problem has also found applications in several 
other areas where pairing correlations are stronger than 
in the BCS theory, e.g., the case of a strongly interacting 
Fermi gas with large scattering length. Let us consider 
the grand canonical Hamiltonian, K ~ H ~ fxN, of a 
two-component atomic Fermi gas, 

k,<T 

UkW t t /-IN 

+ ^j7~"k+iqT"-k+iqi"-k' + iqi"k'-HiqT.UJ 

q,k,k' 

where aj^.^ , ako- are the creation and annihilation oper- 
ators for free atoms with momentum Kk, energy ek = 
?i^fc^/2m, the magnitude of the momentum k = |k|, with 
the pseudospin index a =1,1 denoting one of the two 
trapped hyperfine spin states, and with m and /i being 
the mass and chemical potential of the atoms, respec- 
tively. The single atom kinetic energy measured with 
respect to the chemical potential is Ck = £k — ^J', and the 
momentum dependent Fourier transform of the real space 
inter-atomic interaction is J7kk' = U{\'k — k'|) which in 
the crossover model with a bare attractive contact inter- 
action has the constant strength, [/kk' = U. Due to the 
singular nature of the contact interaction unphysical di- 
vergences arise in three dimensions. The divergences may 
be removed by a standard regularization procedure |2^ . 
By relating the bare potential to the two-body scatter- 
ing length ap and the corresponding effective strength, 
= 4:Trh^aF/iTi, the bare potential can be eliminated in 
favor of the scattering length. This is formally accom- 
plished by using the relation between bare and renormal- 
ized strengths, = — ^J2\k\<k (2ek)~^, where 
V is the volume of the system, and hkc is a momen- 
tum cutoff. At fixed of it is a relation between U and 
kc. The final results become finite and fcc independent 
by taking the limit kc ^ oo aX the end of the calcula- 
tion. The l/£k term regulates the divergence of the free 
two-particle propagator. A similar procedure has been 
applied to the general case of model potentials with non- 
zero range and is presented in detail in the next section. 
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Figure 1: Diagrammatic structure of the equations for the 
dressed finite temperature propagator G (thick arrow lines). 
First line is the Dyson equation for the atomic propagator, the 
second line is the self-energy (shaded disc) derived from the T- 
matrix approximation, where the T-matrix (shaded square) is 
found from the integral equation in the third line. To (white 
square) is the effective two-body interaction which is found 
after proper renormalization of the bare two-body interaction 
U (wiggly lines) by inverting the equation for To in the fourth 
line. The non-interacting propagators are represented by the 
thin arrow lines, and the internal pair of fermions in the third 
and fourth lines are the regularized pair propagator, Xt find 
non-interacting pair propagator, xo, respectively. 

A. Nonzero temperature Dyson equations 

In the standard crossover model in E E Ei the 
BCS gap equation is supplemented by an equation for the 
number density. At fixed particle density this is an equa- 
tion for the chemical potential, which in the BCS limit is 
the Fermi energy, n « Ep. In the BEC limit it is half of 
the bound state energy, fi = —Eo/2, with Eq = /ma^p. 
In the nonzero temperature Green's function formulation 
El of the crossover model the appearance of the bound 
states with finite center-of-mass momentum, q, at suffi- 
ciently large interaction strength is accounted for by use 
of the T-matrix approximation for the self-energy. The 
T-matrix approximation diagrammatically corresponds 
to a resummation of an infinite series of particle-particle 
ladder diagrams. By using methods of nonzero tempera- 
ture Green's function theory the equations for the tran- 
sition temperature and critical chemical potential are de- 
rived from the corresponding Dyson equation within the 
T-matrix approximation, 

G(k,c^„) = Go(k,c^„)+Go(k,w„)I](k,c^„)G(k,c^„), (2) 

with k being an atomic wave- vector, w„ = (2n + 1)ti/ (3 
is the fermionic Matsubara frequency for integer n, and 
where Go and G are the noninteracting and dressed finite 
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temperature Green's functions for the atoms, respectively 
(see Figure m for a diagrammatic representation). Unless 
specified we use below Ti — 1, ks — 1, and (3 = l/{kBT) 
for the inverse temperature. The self-energy S is eval- 
uated in the ladder approximation and is conveniently 
expressed in terms of the T-matrix which is the sum of 
ladder diagrams. In compact notation the equation for 
the many-body T-matrix T is T = U + U xT, where U 
is the momentum dependent bare interaction potential 
and X is the in-medium pair propagator. The product 
UxT is to be understood as a momenta convolutions. 
The convolutions will in general give rise to angular in- 
tegrals. However, in relation to Feshbach resonances the 
important physics resides in the energy dependence of 
the T-matrix to which we make a low energy expansion. 
The angular integral can therefore be applied to the pair 
propagator alone and the angular dependence vanishes 
from the T-matrix equation. In general the T-matrix de- 
pends on the two incoming particle momenta, ki, k2 and 
two outgoing k']^ , kg . For the off-shell T-matrix there are 
in addition dependence on two incoming uji , 102 and two 
outgoing and ^^,^2 particle energies. It is convenient 
to pass to relative momenta k, k' and energies w, u)' and 
total momentum and energy, q, Q.. For the on-shell to- 
tal energy, SI = wi -f = Efei + ^fca = 2£a: + ef , with 
the free fermion dispersion relation, eu — fi^k"^ / (2m), the 
relative momentum in center-of-mass frame is fc = |k|, 
the center-of-mass momentum is g = |q| and the free 
pair dispersion is = ?i^g^/(2mB)with tub = 2m. The 



analytically continued, iVlu 
propagator is 



f2 -|- iO^, in- medium pair 



X(k,q,r!) = i^Go(e+,f^)Go(e-,f^- 



(3) 



where ^± = £±k+q/2 — A*- We tacitly used the fact that 
the T-matrix is independent of the relative Matsubara 
frequencies. Let the two-body limit of the T-matrix be 
denoted by Tq and the free two-particle propagator by 
Xo- The corresponding two-body limit of the many-body 
T-matrix equation then becomes Tq = U + UxqTo- The 
renormalization procedure consists of re-expressing the 
T-matrix in terms of two-body quantities which on the 
other hand can be expressed directly in terms of the scat- 
tering phase shift of the two-particle wave function in 
relative coordinates. Heuristically, we may define L from 
the equation U = {1 — xo)To = LTq and therefore find 
T = L-^LT = L-\T - xoT) = L-^U{1 + (x - Xo)T). 
Hence, after eliminating the bare interaction the T- 
matrix equation becomes 



for small momenta the function g{k) = fccot^o can be 
expanded in relative energy ~ h^k'^/m. Thus, g{k) = 
—ap^ + rck^/2 + . . . , where ap is the scattering length 
of the fermions, and reis the effective range. The shape- 
independent approximation applied here consists of us- 
ing only the two first expansion coefficients which char- 
acterizes the potential by its depth and range, but is 
more generally also a valid approximation for multichan- 
nel scattering problems. For model potentials the shape 
dependence is reflected in the signs and magnitudes of the 
higher order expansion coefficients. The effective range 
has dimension of length and is for an attractive mono- 
tonic potential at resonance, ap ±00, the range of 
the inter-atomic potential. For the resonance phenom- 
ena considered here the effective range can be large and 
negative. Inserting the low energy expansion of g{k) 



into the scattering amplitude, /(fc) = k 



sin((5o) 



{—ap^+rck^/2 — ik) ^, the two-body T-matrix becomes, 
To{k) = -(47raF/m)(l - r^apk'^ /2 - ik)~^ . Adding and 



subtracting 1/26^' in Eq. Q and using the identity 

^^ = ^E(xo(k',q,^^)-2i;). 



(5) 



the free two-particle pair propagator, xoj can be elimi- 
nated from the problem. After inserting the expression 
for ik into the Eq. I^J the inverse effective interaction 
can be defined as 



u:\k) 



1 



1 



-Trk^ 



(6) 



By solving the T-matrix equation the expression for the 
inverse T-matrix becomes 



T-\ci,n)^{u-\E)~x{ci,n)) 



(7) 



where the regularized nonzero temperature in-medium 
pair propagator is. 



X(q,f))-^^(x(k',q,f^) 



1 

2^ 



(8) 



The Matsubara sum in x(k',q, ri) is evaluated by turn- 
ing the sum over Matsubara poles into a contour inte- 
gral using standard techniques [l^l . The regularized pair 
propagator x then becomes 



k' 



nFi^+)-nFii--n) 1 



2£k 



■(9) 



T = To + To{x-Xo)T. 



(4) 



The off-shell two-body T-matrix To(fc) is related to 
the two-particle off-shell scattering amplitude f{k), as 
Toik) = —Airfi^ f{k)/m — —[g{k) + ik]~'^ , where we 
introduced the relative momentum of the two collid- 



ing atoms, k 



and where 



where npiO = i^^^ + is the Fermi-Dirac distribu- 
tion function. At a later stage we are going to make 
a low q, n expansion of y and therefore keep the fl in 
nF in the numerator [2J|. For certain values of (3 and 
/i the pair propagator x can be calculated analytically, 
but it is generally evaluated by numerical integration. 
The angular dependence x = cos0 = k' • q/k'q is absent 
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in the denominator and the angular integration can be 
carried out analytically in the numerator. The regular- 
ized pair propagator thus reduces to a one dimensional 
integral over the magnitude of the relative momentum. 
After renormalization and making the effective range ex- 
pansion the equation for the T-matrix is of the same form 
as in the corresponding equation obtained from the bare 
contact potential with the only difference being the en- 
ergy dependence of the effective interaction, C/* U^,{E) 
where E — h'^k'^/m. Having obtained the expression for 
the T-matrix the self-energy becomes 

E(k, u;„) = ^J2 T^(q. ^^.)Go(q - k, 17, - c^„), (10) 

where ili, — 27rn//3 for integer i' are the bosonic Matsub- 
ara frequencies. 

To locate the transition temperature we apply the 
Thouless criterion j23] which states that the transition 
to the superfluid state is found from the singular be- 
havior of the T-matrix in the long wavelength and low 
energy Hmit. At Tc the Thouless criterion is equiva- 
lent to the gap equation. For a fixed number of atoms 
and at the transition temperature the crossover equa- 
tions are then T~^{q,rt) = for q ^ and ^ to- 
gether with the equation for the particle number conser- 
vation. In order to calculate the particle number density 
n = N/V self-consistently we use the thermodynamic re- 
lation, N = —df^il, where ft is the grand canonical poten- 
tial, e~^^ = Tr [e~^^] , tracing over all degrees of free- 
dom. At the transition the result for the normal phase is 
still valid and the thermodynamic potential can be found 
from the Luttinger-Ward expression for the grand poten- 
tial 

n[G,n = m-^ E e^^^^°"[ln(-Goi-|-S)-fSG], 

(11) 

where the physical dressed Green's function G(k, uj) for 
the atoms and the self-energy S(k, w„) are found from 
the stationary points of ft with respect to variations of 
both G and E. Here, Sfl/ST, = yields the Dyson equa- 
tion Eq. ^ and 5fl/5G = links the self-energy, S, to 
the functional <& which is the generating functional for 
the skeleton graphs and is here constructed in order to 
generate the ladder approximation by functional differen- 
tiation, E — ^S^/SG [2^. In the ladder approximation 

A self-consistent treatment of the zero effective range 
crossover problem has been presented in [21, 21,^]. We 
have employed a non-selfconsistent approach where the 
dressed Green's functions in the pair propagator and the 
self-energy has been replaced by Gq and which is equiva- 
lent to the functional integral approach of 0|. In this ap- 
proximation the grand canonical potential takes the form. 



n = n^+Sn, with = Ek,..„,. e*"''°^ lnGo(k, ujn)/f3V 
and Sn = J2c^,n^ e^''"°^ lnT(q, n,)/pV 

III. TRANSITION TEMPERATURE 

In a BCS superfluid the dispersion relation for the 
quasiparticle excitations is, i5k = (Ck + ^^)^^^' which 
displays an energy gap, A ~ U{tp'4')j with (V'V') being a 
anomalous thermal correlator. Due to the energy gap in 
the single particle spectrum thermodynamic and trans- 
port quantities, such as the heat capacity and thermal 
conductivity, may become exponentially suppressed in 
the superfluid phase and may therefore be used as ther- 
modynamic signatures of the phase transition. Apart 
from such quantities being difficult to observe directly in 
a trapped atom gas there are in a strongly interacting 
Fermi gas additional contributions to the quasiparticle 
gap which may further complicate the identification of 
the superfiuid transition from such observables. A sim- 
ple way to see the appearance of such a non-BCS gap 
contribution at zero temperature is to note that the gap 
in the quasiparticle spectrum is the minimum of i?k over 
all k's. In the BCS limit the chemical potential is positive 
and approximately the Fermi energy, and the gap is A 
for k — kp. By contrast when the chemical potential be- 
comes negative on the BEC side the minimum is shifted 
to fc = for the gap {y? -\- A^)^/^ Q. The appearance 
of /X in is related to the molecular degrees of freedom 
which also contribute to the gap at higher temperatures. 



A. Scattering length only 

Let us initially consider the zero range limit, = 0, 
which results in the standard BCS to BEC crossover of 
a Fermi gas with contact interaction |15, Ji| . The tran- 
sition temperature Tc — l/Pc and the chemical poten- 
tial fic are at the transition determined from the Thou- 
less criterion, \imq^Q^oT~^{ci,V,) = 0, together with the 
number density equation, n = jj.) +6n{f3, fi), where 

the fermion contribution is n^(/3,/^) = J2k <7^Fi^k)/^^- 
The fiuctuation term Sn derived from the ladder self- 
energy contribution to the thermodynamic potential is 
5n(/3,M) = -EqEa„5Mlnr(q,0,)//3V. To simplify 
the calculation of the bosonic number density we follow 
^EI^HMHEI^ and make a low energy and long wave- 
length expansion of the inverse T-matrix 

T-\q,n) = zin-e^-^iB), (i3) 

where ifli, — > fl before frequency expansion, the molec- 
ular dispersion relation is = ?i^q^/2m|j with rrig 
being the renormalized boson mass, the wave function 
renormalization constant Z = Z{(3,^), the effective bo- 
son mass m|j = mg{P,fj,) and the bosonic chemical po- 
tential idB = A'b(/3jM) are found from the low order en- 
ergy and momentum expansion coefficients of x together 
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Figure 2: The transition temperature for re = on the BEC 
side of the crossover is plotted as a function of the inverse 
coupling constant from the full calculation (full curve) and 
the corresponding result (dashed curve) obtained from the ex- 
pression for the transition temperature for the non-interacting 
BEC with the numerically calculated effective boson mass. 
The inset shows the crossover behavior of the effective bo- 
son mass in units of the fermion mass in the whole crossover 
region. 



with U^r^{k). Re-expressing the frequency integration as 
a standard sum over bosonic Matsubara frequencies Sn 
becomes 



(14) 



I) 



is the 



where = - ^b, and nB{£,) = 
Bose-Einstein distribution function. 

Before we turn to the results of the numerical cal- 
culations we note that the results for the uniform case 
with a contact potential are well-known: In the weak 
coupling limit as (fc^aF)"^ —oo the Fermi gas is ex- 
pected to make a transition at the transition temperature 
Tc = (8e'^-V7r)ri.e-'^/2feF|a^| ^ pair-correlated BCS- 

state, where the Fermi temperature is, Tp = EpjkB, 
and 7 — 0.57721 ... is the Euler-Mascheroni constant. 
For a uniform noninteracting two-component Fermi gas 
with the Fermi energy Ep the corresponding Fermi 
wave vector is kp = {2mEp /fi^Y/"^ , and the parti- 
cle density is n = kp/Sn'^. As the inverse scatter- 
ing length is increased strong coupling effects appear 
and when [kpap)'^ passes through zero from below 
a two-body bound (molecular) state appears, and the 
fermionic states hybridize with the molecular states. 
When the inverse scattering length is increased further 
{kpap)~^ — > oo the Fermi surface smears out (which 
can be seen by inspection of the momentum distribu- 



tion) and a molecular Bose gas is formed which under- 
goes Bose-Einstein condensation at the transition tem- 
perature Tbec = (?iV2mfcs)27rC(3/2)-2/3(^/2)2/3 ^ 

0.2181^. Between the two limits the system smoothly 
crosses from BCS-like behavior to BEC-like behavior in 
the region, —1 < (kpap)^^ < 1. In the BEC limit where 
the fermion contribution vanishes Sn reduces to twice 
the boson density obtained from the Bose-Einstein dis- 
tribution, Sn{f3,fi) = 2J2qnB{tg)/V, for bosons with the 
mass twice that of the atoms, m*g — tub = 2m. 

The problem of numerically solving the pair of equa- 
tions n = riQ + Sn and T~^(0,0) = for f3c and fic was 
solved iteratively using multidimensional root solvers im- 
plemented in [^. At each iteration the momentum inte- 
grals appearing in to|j, Sn, ^b and in the fi derivatives 
appearing in the number equation were numerically inte- 
grated using appropriate adaptive quadrature algorithms 
[231 also reimplemented in |32 |. The numerical results for 
the zero effective range crossover are presented in Fig. 
(full curve) and show a smooth crossover with charac- 
teristic Tc peak near the zero of the chemical potential 
and qualitatively in agreement with previous calculations 

For a weakly interacting Bose gas there are a num- 
ber of corrections due to critical fluctuations close to Tc 
with the leading term being linear and positive in the 
bosonic gas parameter cx kpap and indicating a peak 
in Tc- Mean field corrections to the transition temper- 
ature for a weakly interacting Bose gas with nonzero 
effective range give rise to a leading correction term 
which is negative and cubic in the gas parameter and 
thus indicates absence of peak in Tc. However, the non- 
selfconsistent Gaussian model presented here includes 
none of the above mentioned boson-boson effects and the 
correction to Tc observed on the BEC side here is due 
to an externally induced change of the molecular binding 
energy Eq. When {kpap)~^ is lowered by increasing the 
strength of the effective interaction the binding energy 
goes down and the molecules dissociate more easily at a 
fixed temperature and in addition the fermion states be- 
come populated. The coupling to the fermions also gives 
rise to a renormalized boson mass, < niB, where the 
bare boson mass is = 2m. Therefore the ideal Bose 
gas condensation temperature is shifted to higher values, 
Tc = Tbec + STbec where STbec/Tbec cx \SmB\/mB 
for decreasing effective mass, SmB — nn*^ — tob as 
{kpap)~^ increases. In Figure [2 the transition temper- 
ature in the BEC Hmit is plotted as a function of the 
{kpap)~^ using the numerically calculated Tc (full curve) 
and using the ideal Bose gas expression for the transition 
temperature with T* = T'^p.(j oc 1/ m*^ using the numer- 
ically calculated effective mass m*^. The inset in Figure[2l 
shows the effective mass m*^ identified from the asymp- 
totic behavior of the T-matrix, see Eq. . The effec- 
tive mass becomes as expected twice the fermion mass in 
the BEC Hmit, but no attempt was made to decompose 
Sn into partial contributions. Thus, the effective boson 
mass becomes a hybrid object in the crossover region 
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Figure 3: The BCS-BEC crossover transition temperature of 
the homogeneous Fermi gas as a function of the inverse cou- 
pling constant, —2 < (kpap)'^ < 2, for several values of the 
effective range parameter kprc- The zero range result (full 
curve) corresponds to the usual Tc curve for the crossover. For 
increasing values of the —kpre the transition temperature is 
suppressed in the crossover region. In the inset we character- 
ize the small effective range to large effective range crossover 
by plotting the half-density crossing point in {kpap)'^ as a 
function of the —kpre- The curve shows a shift in the position 
of the crossing to lower couplings for increasing —re. 

without clear physical interpretation as a single particle 
mass. 



B. Effective range expansion 

As indicated in the introduction narrow Feshbach reso- 
nances can be mimicked by an interaction potential with 
large and negative effective range. To account for such a 
situation the zero effective range contact potential must 
be replaced by a potential with nonzero range. How- 
ever, for monotonic and attractive potentials the effective 
range is at resonance the averaged microscopic range of 
the potential and thus positive [2J|- For potential models 
the effective range can be related to the structure of the 
interaction potential and can be negative in the presence 
of repulsive cores or barriers. In the latter case consid- 
ered here the negative effective range is related to the res- 
onance physics of potentials with barriers. For a poten- 
tial with a barrier quasi-bound states may appear which 
have a finite lifetime and thus decay into the continuum 
states. In the zero temperature calculation [rRl3fl| it was 
suggested to mimic the presence Feshbach resonances by 
using a two-body square well with a square barrier model 
potential in place of the contact potential. In contrast 
we use the renormalized crossover equations of Section 



mi and the numerics of subsection IIII Al to examine the 
effect of the energy dependent effective interaction 

U^k) = -^^fl-lr^apkA'' (15) 
m \ 2 J 

on the transition temperature for values of the effective 
range —kprc ranging from small to large and positive. 

In Figure the numerical results for the tran- 
sition temperature for several values of —kpr^ — 
1,2,4,8,16,32,128 are plotted as a function of the 
{kpap)~^ together with the zero range = result 
(full curve). For a broad resonance with vanishing kpr^ 
the usual zero effective range result is obtained. In the 
crossover region the transition temperature is sensitive 
to the effective range and the peak in Tc becomes sup- 
pressed for increasing values of —kpTc- In the BCS (BEC) 
limit a negative shifts the effective interaction to lower 
(higher) values corresponding to a lower transition tem- 
perature which is consistent with the decrease in Tc for 
increasing —r^ at fixed ap- Furthermore, it is seen that 
the energy dependence of the effective interaction is only 
important in the crossover region for —kpVc ^ 1. The 
small value of the parameter (fci?ai?)(fcF''c) the energy 
dependence of C/* is suppressed and therefore both the 
asymptotic BCS and BEC limits are recovered even for 
large kprc- 

It is tempting to compare our result to the corre- 
sponding calculation of Tc for the atom-molecule Hamil- 
tonian. In j^H H3 the transition temperature was cal- 
culated as a function of the detuning v for a homoge- 
neous Fermi gas both for weak inter-channel coupHng 
9bf — O.GEp and for strong coupling gsF — 20Ep. 
In addition to an explicit contribution to the number 
equation arising from the bosonic degrees of freedom, 
the effective interaction appearing in the Thouless cri- 
terion also includes the effect of repeated fermion-boson 
scattering in addition to the fermion-fermion scatter- 
ing, U*{q,n^) = Uff - g%pDo{q,n„), where Uff is 
the inter-fermion coupling strength, gsF is the boson- 
fermion coupling strength, ZJ^f ^(q, fii/) — (iil^ — ^^)^^ is 
the free molecule propagator, = h'^q'^/ {2mB) — 2{^—h') 
is the bare molecule dispersion relation. For small cou- 
plings corresponding to a large and negative effective 
range Tc was found to increase monotonically from the 
BCS to the BEC limit. In the opposite limit of large 
couplings corresponding to a small effective range a pro- 
nounced peak close to threshold appears. These results 
are consistent with the trends found in the present calcu- 
lation. It is also interesting to compare the result to the 
transition temperature obtained for a mixture of fermions 
and long-lived bosons (narrow resonance case) with the 
total density constrained by total fermion number con- 
servation, and with tob = 2m. By decomposing the total 
density into three contributions, n = Uq + 2n^ + 6n, and 
by neglecting 6h when solving the crossover equations, a 
monotonically increasing Tc is found which appears to be 
the Tc curve in the large —r^ limit. The effective range 
induce a cutoff in the momentum integrations and thus 
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reduces the influence of high momentum modes. This 
decreases the effect of the fluctuations and leads to a 
suppression of Tc which is again consistent with the in- 
terpretation of the large Tc peak being due to Gaussian 
fluctutations. To characterize the 'crossover' from small 
Te behavior to large behavior we notice that at a flxed 
kpap the fermion contribution to the total density be- 
comes suppressed for increasing —r,, due to decrease in 
the chemical potential in the crossover region. Therefore, 
the half-density crossing point of the partial densities, 
= (5n = 5, is shifted to the BCS side. The inset in 
Fig. El shows a semilogarithmic plot of the half-density 
crossing point as a function of kpTc- As value of the ef- 
fective range is increased the crossing point shifts from 
{kFaF)~^ ~ —0.07 for rc = to {kpap)'^ ~ —0.5 in the 
large —Tc limit where the dependence on the effective 
range is expected to drop out. Therefore, the transition 
temperature appear to have a lower bound at flxed cou- 
pling for large —r^ when neglecting high order terms in 
the expansion of k cot S. 

It is to be noted that the current calculation does not 
account for i) the induced interactions due to vertex cor- 
rections which reduce Tc by a factor of 2.2 in the BCS 
limit ^|. ii) the residual boson-boson and boson-fermion 
interactions which are characterized by the boson-boson 
scattering length, as, and the boson-fermion scattering 
length, ubf, on the EEC side as recently found from 
solving three and four-body problems [2^. Let us flnally 
mention that to the extent that many-body physics is 
important to the crossover problem it must be stressed 
that many theoretical issues, e.g., non-zero effective range 
calculations, renormalization of coupling constants, in- 
clusion of higher order diagrammatic contributions, and 



generalization to pairing states in higher angular momen- 
tum channels, are easier to handle within the single chan- 
nel potential model. 



IV. SUMMARY 

We presented the renormalized T-matrix equation 
which allowed the non-zero effective range BCS-BEC 
crossover equations at the transition to be formulated 
directly in terms of two-body scattering properties, the 
scattering length and the effective range. By matching 
the scattering properties of the generalized single-channel 
BCS-BEC crossover model with the two-channel Fesh- 
bach resonance atom-molecule model it was argued that 
the internal energy scale of a resonance translates into to 
the effective range of the potential model. Broad reso- 
nances correspond to a small effective range and for nar- 
row resonances the effective range become large and neg- 
ative. By using the NS crossover equations with the self- 
energy derived from the ladder approximation we numer- 
ically calculated the superfluid transition temperature of 
the generalized BCS-BEC crossover model as a function 
of the two-body scattering length and found Tc to be 
suppressed in the crossover region for increasing values 

of —kpTc- 
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